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1. Introduction 



One of the first observables studied in QCD was the energy-energy correlation (EEC) 
in e + e~ annihilation [1, 2, 3]. It is a function of 9, the polar angle between pairs of 
hadrons weighted by their energy. A similar correlation in DIS involves not a polar 
but an azimuthal angle. Here the observable is given by 



H(x) =sinx^ Pth ^2 h ' 5(cosx+cos(f) hh ,) , 
hh' V 



»hh> 



= <j>h-<t>h>, o<x<tt, (i.i) 



where the sum runs over all outgoing hadron pairs with p t h and <ph the transverse 
momentum and azimuthal angle in the Breit frame. We take — n < tp^y < ir. The 
distribution for the observable H(x) is dominated by events with large p t jets: indeed 
we wish to exclude events with p t much less than Q in order that we can perform the 
calculation based on a dijet type configuration. On the other hand, to simplify the 
QCD analysis and consider the incoming 7 or Z as pointlike, one should limit the 
distribution to events in which all jet p t are not much larger than Q. A convenient 1 
way to set these conditions is obtained by limiting the dijet resolution variable 7/2 
defined by the k t jet finding algorithm [4] as follows: 

y- < y2 < y+ ■ (1.2) 

Here the lower bound y_ selects large p t events, while the upper bound y + allows us 
to consider the incoming vector boson as pointlike. 

The azimuthal differential distribution for this observable is then defined, at fixed 
xb and Q 2 , by 

da{x,y±) _ ST f da n 



dx dx-^dQ' 



with dcr n /dxBdQ 2 the distribution for n emitted hadrons in the process under con- 
sideration. One introduces the following normalized azimuthal correlation (EtEtC) 

d£(x,y±) = -1/ \ da(x,y±) 
d X a {y±> d X dx B dQ 2 ' 

-(^) = E / dx^ 2 eiy+ - m)e{m ~ y - ] f ^ 

By taking into account the case h — h! in (1.1), i.e. x — n > this correlation is 
normalized to 1. We will be concerned with the back-to-back region of small x — 

71 - \(f)\. 

1 This method to set a limit on the large transverse momenta should contain small hadronization 
corrections, see [5, 6, 7]. 
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A standard method to start to understand the behaviour of this correlation for 
small x consists in trying to exponentiate the leading logarithmic one-loop contribu- 
tion given 2 by 

n> 1 

£(x,y±) = l-C r /m 2 - + ..., C T = 2C F +C A , (1-5) 

where Ct is the total colour charge of the three hard partons qqg generating both 
the outgoing dijet and the incoming proton jet. The exponentiation of this one-loop 
term gives a characteristic Sudakov form factor leading to an azimuthal differential 
distribution dY*/dx with a peak at small x near 

X ~ e -T/(2C T as) _ ( L6 ) 

This is actually what happens for EEC in e + e _ annihilation in the back-to-back 
region (small X = n ~ \@\ with 9 the polar angle): the leading logarithmic one-loop 
contribution of S e + e -(%), is given by the same expression in (1.5) with Ct = 2Cp, 
the charge of the qq pair in the hard vertex. In this case the complete QCD analysis 
[1, 2, 3] and the data [8] confirm the presence in dT> e + e -/dx of the Sudakov peak in 
(1.6). 

There are experimental indications that azimuthal distributions behave differ- 
ently. For instance, in pp processes, the differential distribution in the azimuthal 
angle between two hadrons behaves smoothly in the back-to-back region and shows 
no sign of a peak, as reported for instance in [9]. We should then understand why 
EtEtC and EEC behave differently. 

In this paper we perform the QCD analysis of the DIS azimuthal correlation 
(1.4) in the back-to-back region. The analysis is similar to the one [1, 2, 3] for EEC 
with a few major differences: 

• the distribution (1.4) is dominated by dijet events which originate from the 
qqg hard system (3- jet events including the beam). Only recently has the QCD 
analysis been extended to distributions involving 3 jets [6, 7, 10] at the same 
accuracy as available in the study of 2-jet CIS distributions [11, 12]; 

• the radiation off the incoming parton contributes not only to the observable 
H(x), but also to the parton density evolution; 

• QCD resummation is performed in the 1-dimensional impact parameter b con- 
jugate to the azimuthal angle. In the EEC case the variable is a polar angle 
and then one deals with the standard 2-dimensional impact parameter. 

As in other similar analyses, see [7, 10, 12], due to QCD coherence, the radiation 
factorizes giving rise to the evolved standard parton density functions times the "ra- 
diation factor", a collinear and infrared safe (CIS) quantity involving 3-jet emission; 

2 Contributions from parton density functions in E T EiC enter beyond double logarithmic level. 
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We work at the following accuracy: (i) double (DL) and single (SL) logarith- 
mic perturbative (PT) resummation (for the logarithm of the distribution in the 
^-representation we resum all terms a™hi n+1 bQ and a™ln n bQ); (ii) matching of re- 
summed with exact fixed order results (in this paper only to one- loop order); (iii) 
leading non-perturbative (NP) corrections coming from the fact that the QCD cou- 
pling runs into the infrared region [5, 13]. 

The paper is organized as follows. In section 2 we discuss the QCD process and 
the observable H(x). In section 3 we describe the procedure for the resummation 
at SL level. In section 4 we present the result of the resummation and evaluate 
the leading power correction. In section 5 we present the numerical evaluation. In 
section 6 we the summarize and discuss our result. Appendices contain the technical 
details. 



2. QCD description 

We work in the Breit frame of the DIS process with P and q the momenta of the 
incoming proton and the exchanged vector boson (7 or Z ) 

, = §(0,0,0,2), P= ^(1,0, 0,-1), *b = ^. (2.1) 

We first discuss the QCD hard elementary vertex and then the secondary parton 
emission. 

2.1 Elementary hard process 

DIS dijet events originate from the elementary hard vertex 

qPl ^P 2 P 3j (2.2) 

with Pi along the z-axis. We take P 2 , P3 in the {y, z}-plane, the "event-plane", see 
appendix A. We use the kinematical variables 

{= w x= w^j >XB - (2 ' 3) 

The outgoing transverse momentum P t and the squared invariant masses Q 2 ab = 
2(P a P b ) are given by 



Pt=QJt(i-Z) — , Q 2 n = -Q\ QL— Q 2 , Qls = —Q 2 - ( 2 - 4 ) 

V x x x x 

We distinguish P 2 from P 3 by assuming 

(P X P 2 ) < (PxP 3 ) , (2.5) 
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which restricts us to the region < £ < \. In appendix A we report the elementary 
distributions for the process (2.2) as functions £ and x. They depend on the "con- 
figurations" of the three partons. If we consider only photon exchange it is enough 
to identify the gluon momentum, we use the index 5 = 1, 2 or 3 to denote that the 
gluon momentum is P\,Pi or P 3 respectively. If we consider also Z exchange, with 
parity-violating terms, we need to identify also the nature of incoming parton, this 
will be identified by the index r = q, q,g. Finally for quarks we need to introduce 
also the flavour index /. The configuration of the hard vertex will be then defined 
by the configuration index p = {5, r, /}, see appendix A for a detailed presentation 
of these conventions. 

2.2 The observable at parton level 

The QCD process involving multi-parton emission is described by one incoming par- 
ton of momentum pi (inside the proton) and two outgoing hard partons p2,Ps ac- 
companied by an ensemble of secondary partons ki 

qpi ^ p 2 P3h- • -K , p 1 = x 1 P. (2.6) 

Here, taking a small subtraction scale \x (smaller than any other scale in the problem), 
we have assumed that p\ (and the spectators) are parallel to the incoming proton. 

Since H(x) in (1.1) is linear in p h and py we may replace the sum over hadrons 
with a sum over outgoing partons for the process (2.6). To compute the distribution 
for small x & t SL accuracy, we can consider the secondary partons to be soft and 
take the contributions to H(x) up to the first linear correction in the soft limit. Hard 
collinear radiation in the jets gets embodied into the hard scale of the corresponding 
radiator, see [11], and, in the case of the incoming jet, contributes to the parton 
density function, see [7, 10, 12]. In this limit the outgoing hard partons p 2 and p 3 
tend to P2 and P 3 , with soft recoil we need to consider to first order. From appendix 
B, we have 

W^^ll-^j, X 2 3 = vr-|0 23 |^|0,|, ^ = E^' (2 ' 7) 

with k xi the out-of-event-plane momentum component. The event plane is defined 
in appendix B and P t in (2.4). Neglecting contributions to H(x) from two secondary 
soft partons we have (see figure 1) 

^(x)-^|(^1-E^)^-^3) + E|[^-^) + ^-^)]}- (2-8) 
From the expression of Xai = n — \<t>ai\ we obtain in the soft limit (see appendix 

B) 

[Kx-X2i)+5{x-X?>i)\ - 5 (x - \4>i - 0x|) , \k V i\ - ^ I cos 0,| , (2.9) 
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Figure 1: The kinematics of azimuthal correlation in the plane orthogonal to the Breit 
axis. The hard partons p 2 and p 3 correspond to the configuration 5=1 and ki is a soft 
secondary gluon. 

with (pi the azimuthal angle between ki and either the hard parton p 2 or p 3 : 

(f) i2 for ki near p 2 , 
3i for ki near p 3 . 

Splitting the expression of H(x) m to a 'hard' and a 'soft' contribution we can write, 
to first order in softness, 



H( X ) ^ H h ( X ) + H S ( X ) , (2.11) 

where, 

IP 2 

Mx) = -^Kx-\<i> x \), 

2P 2 ^ fo. ( 2 - 12 ) 

This form explicitly shows that H(x) is, by itself, a CIS observable in the sense 
that the azimuthal correlation has nothing more than the singularities of the total 
DIS cross section, namely the incoming parton collinear singularities leading to the 
anomalous dimension, see later. Indeed, the hard term (hard-hard contribution 
with no recoil, pt2Pt3 Pf) depends only on the total recoil \<fi x \ from secondary 
partons. The soft term H s (hard-soft contribution together with the recoil piece 
of hard-hard contribution) is present only in the real contribution. Here the soft 
singularity of the matrix element is damped by the ku factor, while the collinear 
singularity, 0«^O, is regularised by the vanishing difference of the delta functions in 
the square brackets. 
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3. Resummation 



Considering the region x < 1, the starting point for the QCD resummation of the 
azimuthal distribution is the factorization [14] of the square amplitude for the process 
(2.6) with n secondary partons soft or collinear to the primary partons P a in the hard 
vertex (2.2). For each configuration p we have 

\M n (h ...k n )\ 2 ~ |M | 2 • S n (h ...k n ). (3.1) 

The first factor is the squared amplitude for the elementary vertex (2.2) for the 
different configurations p. It gives rise to the elementary hard distribution da p as 
function of the kinematical variables x and £ in (2.3), see appendix A. The second 
factor is the distribution of secondary partons emitted from the hard momenta P a . 
Since in the emission, the nature of incoming parton changes, S n is a matrix in the 
configuration index p. Corrections to the factorized expression are relevant for finite 
X and can be taken into account by a non-logarithmic coefficient function. From 
(3.1) the azimuthal correlation (1.3) is then given by 

da{Xi V±) f XM dx [t+ ( db p 



dx dxB dQ 2 



where the limits %, £± are determined by the values of y± which limit the transverse 
momenta, see appendix A. Here C p (a s ) is the non- logarithmic coefficient function 

Cp (a s ) = l + ^c[ p \ X ) + 0(a 2 s ) . (3.3) 

71 

The distribution I depends on the geometry of the dijet system through the hard 
scales Q ab given in (2.4). It is given by 

I( X )=J Q dx 1 V(x ul i) J dr n -S n (h...k n ), (3.4) 

where V{x\,p) are the parton density functions at the (small) subtraction scale p and 
x\ is the incoming parton momentum fraction (2.6). The phase space and observable 
are contained in 



dT n = TT / ^ • {H h ( X ) + H S ( X )} •*(—-*! TT *) ■ (3.5) 

t=i J nuJi \ x hi ) 



The last delta function accounts for the fact that the incoming parton pi loses en- 
ergy by radiation. Here Zi is the collinear splitting energy fractions for radiation of 
secondary partons fcj in the region C\ collinear to p\. Note that it is this last delta 
function which is responsible for the non cancellation of the (residual) incoming 
parton collinear singularities which build up the parton density function. 
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The resummation of the distribution X is based on the factorization of soft and 
collinear contributions in S n . Then one needs to factorize also the expression of dT n . 
For the observable H(x) this is done by 1-dimensional Fourier transform and, for the 
delta function in Xb, by Mellin transform 




where the last integral runs parallel to the imaginary axis with Re(iV) > 1 and 

e( Zl ) = zf~\ for ki e d , e(z t ) = 1, for h £ & . (3.7) 

Each secondary parton contributes to dT n by a factor 

[6(^)008(6^)-!] = [e(zi)-l]cos(bk xi ) - [l-cos^/c^)] , (3.8) 

with the one included to take into account virtual corrections. Since S n is even in 
all k X i, we have replaced e lbkxi cos(bk X i). 

The splitting made in (3.8) shows that we can resum independently the parton 
density function and the radiation function of the observable H. Indeed, from the 
well known SL approximation 

[1-008(6^)] ~ {\k xl \ -r 1 ) , b = e™ \b\ , (3.9) 

we have that the first term in (3.8) involves partons in the region k xi < b^ 1 collinear 
to Pi with weight [z^ -1 — 1]. Their contributions reproduce the anomalous dimension 
and resum to the parton density function at the proper hard scale which is The 
second term in (3.8) involves partons in the complementary region with weight — 1, 
only virtual contributions survive. Their contributions are resummed to give the 
standard radiation factor for the observable H(x) with (virtual) frequencies larger 
than ft -1 . This method to show the factorization of these two contributions has been 
used and discussed in [7, 10, 12]. 

The final result for J thus becomes 

Zt A /(x;Q«6) = t#- / — ™<b x P t ) ■ v Tj (Ab- 1 ) -Mb;Q ab ) , (3.io) 

where the radiation factor is given by the sum of the hard and the soft contribution 
.4,(6) = A h s (b) + A s s (b) , A h s (b) = e- n ^ , ^(6) = ■ B 5 {b) . (3.11) 
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It will be seen that only the hard contribution is required for a PT resummation at 
SL level, while the soft contribution is responsible for leading NP corrections. 
Here TZ$ is the full CIS radiator for the secondary emission and is given by 

/d 3 k 
W s (k) [l-cos(bk x )} . (3.12) 

The quantity B$, representing the H s (x) contribution to the observable, is 

/d 3 k K 
—W s (k) § [cos(6P t 0) - | cos 0|] . (3.13) 

Here Wg(k) is the two-loop distribution for the emission of the soft gluon from the 
primary three hard partons P a in the configuration 5: 

W s (k) = ^ {^Sa + w Sb - j-^Wai}j , a^b^5 , (3.14) 

where w a b(k) is the standard distribution for emission of a soft gluon from the ab- 
dipole. It has a simple expression in the afr-dipole centre of mass. Representing in 
this frame the soft gluon momentum by Sudakov decomposition 

k* = aP: + (3P b * + K, aP=^, (3.15) 



Q: 



ab 



with P* and P b back-to-back momenta {Q 2 ab = 2P*P b *) and k a two-dimensional 
transverse vector, we can write 

a s (n) d 3 k d 2 nda k 2 
Wab{k) = =- , = , a > —j- . (3.16) 

UK 1 TXUJ 7T a Q Z ah 

To ensure SL accuracy, we have to take a s defined in the physical scheme [15]. Taking 
P*, P b * along the z-axis, the x-components of the momenta are the same as the ones 
in the laboratory frame with the event plane given by the yz-plane. We then have 

A similar analysis gives the (weighted) total cross section cr(y±) in (1.4) needed 
to obtain the normalized differential distribution dH(x)/dx- We find 

^^£?f<^)fM-.«). 

where the parton density function is reconstructed at the unconstrained DIS hard 
scale Q. 

In conclusion we have that, while in the DIS cross section the proper hard scale 
of the parton density functions is Q, in the distribution (3.10) it is given by b^ 1 . This 
is a result of QCD coherence. 
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4. The result 



The general structure in (3.11) of the radiation factor is similar to the one in EEC: 
the hard term Ag gives the PT contribution at SL accuracy. The soft term Ag is 
beyond SL accuracy. On the contrary, for the NP correction, it is the soft term which 
gives the leading l/Q-power piece; the hard term gives a subleading l/Q 2 piece. We 
discuss these results in the two next subsections. 

4.1 Hard contribution 

We start from the PT part of the radiator (3.12) which is given at SL accuracy by 
the approximation (3.9). This has been computed in [7] (where it corresponds to the 
replacements v — > and vf3 2 , vfiz — > b). We have 

nr(b)=R 5 (b) = J2C^r(b,eQ^) , r(b,Q)= ff^ln|, (4.1) 

o=l Jl /~ b K 71 m 

where the various hard scales and constants are, see [7, 10], 

= Ql 5) = Q ab , cp = cf = C F , tf> = (I s > = e-f , 



n (5) _ QasQsb r (5) _ r AS) _ -M. ( 4 - 2 ) 



with 5 a gluon and a or b a quark (or antiquark), /3 is the first coefficient of the 
beta function (see (D.24)). As is typical for a 3-jet quantity, the scale for the quark 
or antiquark terms of the radiator is the qq invariant mass, while the scale for the 
gluon term is given by the gluon transverse momentum with respect to the qq system. 
The rescaling factors (a '' take into account SL contributions from non-soft secondary 
partons collinear to the primary partons P a . The integration variable k in (4.1) is the 
out-of-plane component of momentum k x and a s is taken in the physical scheme [15]. 
The rescaling factor 2 in the running coupling comes from the integration over the 
in-plane momentum component. The exact expression of hard scales and of rescaling 
factors in (4.2) is relevant at SL level. 

The complete SL resummed result is then given by 

dZ PT ( X ,y ± ) A ^ [*» dx [*+ Jt ( da p 

dx 

where 



17 (x) =[™ "xcotp,) fm • m = 9" v "f 5_i ) e ~ R,m ■ (44) 

The integrand is computed at SL accuracy in 6-space, that is, lnjF p (6) resums all 
terms a™\n n+1 bQ and a"\n n bQ. For x ^ we have that J PT (0) has a l/^/oT s - 
singular behaviour (it is essentially given by the integral of a Sudakov form factor, 
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see [2] for the EEC case). Formally the upper limit of 6-integration is set to infinity, 
but actually, from (4.1) is b < 2/Aqcd- 

As we shall see, the differential azimuthal distribution is monotonically decreas- 
ing by increasing \- This is in contrast with the EEC case in which the differential 
distribution (analogous to dH(x)/dx) develops a peak in the small x region (the 
back-to-back polar angle in this case), which is present also in the data (see [8]). In 
the next section we present the numerical evaluation of the azimuthal distribution 
while for the remainder of this section we discuss the origin of this behaviour and of 
the differences with respect to the EEC case. 

General features To discuss J PT (x) for small x (index 5 neglected here) we divide 
the 6-integral at the point b = l/(xP t ) with X — X elE '- 

A fbo Ah a r°° Ah 

(4.5) 

The sine function oscillates with a first maximum above the splitting point. The 
integrand J-{b) has a Sudakov behaviour at large b (see in appendix C the discussion 
for T in DL approximation). From this behaviour we have that the two contributions 
behave as follows: 

• J( + )(x) is dominated by the lower bound at b = b 

e -RWm) 

I (+) (X) (4.6) 

X 

where, for small x, the remainder S^ + \x) can be expanded in the SL function 

a 1 

R'( X ) = 2C T — In - , C T = 2C F + C A . (4.7) 

7T X 

The Sudakov factor e~ R /x in (4.6) has a peak at x ~ Xi with R'(xi) — 1 which 
corresponds to (1.6). This Sudakov behaviour is expected since large values of 
b force all the emitted k xi to be soft, \k xi \ < xPt- 

• X^~\x) takes contributions from values of b not necessarily large. It is a decreas- 
ing function of x from the value I^~\0) and has a 1/ ^/a^-singular behaviour. 
For b not large, the soft recoil xPt = IX^^I i s obtained by cancellations 
among some larger k xi . 

The relative size of the two contributions in (4.5) can be estimated by the DL ap- 
proximation of J 7 , see appendix C. One finds that the two terms iW(^) become 
comparable at x ~ Xi> i- e - j us t near the Sudakov peak and for smaller x the contribu- 
tion I^~\x) takes over. As a result, the differential azimuthal distribution dE(x)/dx 
continuously decreases from x = and does not have a peak at small x- 
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Comparison with EEC. The behaviour of dS e + e - (x)/dx, the differential EEC 
with x the back-to-back polar angle in e + e~, is very different. This distribution is 
given by an expression similar to (4.4) with the most important difference 3 being 
that in EEC the integration measure is xd 2 b instead of db. As a consequence, the 
contribution involving the finite b regions, which corresponds to X^~\x) in (4-5), is 
relevant 4 only for x^Xi with Xi the position of the Sudakov peak given in (1.6) for 
C T = 2C F . Therefore, the distribution dY> e + e -(x)/dx manifests a Sudakov behaviour 
in the small x region. The fact that the contribution involving the finite 6-region is 
essentially negligible in the EEC case can be argued by observing that it is dominated 
by emitted transverse momenta k t i not necessarily soft which have to sum up to a 
total soft recoil. The constraint of momentum cancellation in 2-dimension is much 
stronger than in 1-dimension as for the azimuthal case. 

A similar way to see the different behaviour in the EEC and E t EtC cases, is 
obtained by trying to factorize the Sudakov form factor e~ R / x from the differential 
distribution. In the EEC case the remainder factor can be expressed in terms of the 
SL variable R'(x) (with C T = 2C F ) in the full range x ~ Xi (with xi the position 
of the Sudakov peak). In the EtE-tC case instead the remainder factor is a function 
of R'(x) only in the range x ~ Xi while at smaller x it does cancel the e~ R /x factor 
and then it is not any more a SL function, see in appendix C the analysis in the DL 
approximation. 

4.2 Full contribution 

To complete the analysis we need to take into account also A s s in (3.11) corresponding 
to H s (x), the soft contribution to the observable, see (2.12). From (3.14) we can write 

B s (b) = ^ (B 5a + B Sb - -^B ab ^j , a + b + 5 , (4.8) 

with the {a6}-dipole contribution given by 

„ f Qab d 2 K aJn) f 1 da k t r „ Tx - n 

B ab {b)= ^Ll I —.± cos&P^ - cos0 . (4.9) 

Jo 71 K l 7T J K 2/Ql b " P% 

Here the integration variables k and a are the momentum components of the soft 
gluon in the {a&}-dipole centre of mass (see (3.15) and (3.16)), while k t and are 
the transverse momentum and azimuthal angle in the Breit frame, see (2.10). In 
appendix D we give their kinematical relations. For instance, for the 12-dipole we 

3 In EEC case the integrand T(b) is simply given by the exponential of a radiator R(b), which 
involves only two hard emitters, the qq pair with total charge Ct = 2Cf, and the cosine is replaced 
by the Bessel function. 

4 From the DL approximation this term becomes comparable to the large &-term, corresponding 
to J(+)( X ) in (4.5), at x - X2 = xl i.e. R'ixi) - 2. See [2]. 
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have 

k y = aP t + K y , tan0 = — - , (4-10) 

Ky 

where we used the fact that the Lorentz transformations between the two frames do 
not involve the x-direction so that k x = k x . The result (4.10) can be easily recovered 
for instance observing that for soft transverse momentum in the dipole frame k — > 0, 
the momentum k becomes collinear to P 2 for a > and to Pi for a — > 0. 

The two terms in the square brackets of (4.9) have DL singularities for k — > 
and 0^0. However the sum (multiplied by k t ) is regular which gives 

A* s (b) ~a 8 - A^(b), (4.11) 

so that this term contributes at PT level beyond SL accuracy, see appendix D and [3] . 
This fact is carefully explained at the end of section 2. Recall that the cancellation is 
between the contributions to the observable H(x) from the recoil of the two outgoing 
hard partons {P2P3} and from a soft secondary and a hard parton {kip 2 /z}. 

The soft part of the radiation factor A s s gives instead the leading NP correction 
which is proportional to b, see also [3] . This NP correction comes from the region in 
which the soft gluon in B ab is collinear to one of the outgoing hard partons P 2 or P 3 . 
For instance for the B 12 radiator, where (4.10) holds, this region corresponds to 

k y ~ aP t > k x , (4.12) 

which gives 

k t ~ aP t , cos(bP t 4>) ~ cos — - , cos ~ 1 . 

a 

Changing variable u = b\n x \/a, and taking in the integrand the leading piece for 
k — ► 0, the contribution to B 12 (b) from this region is 

dtt 2 kch s (k) f w d(f) f°°du f dn 2 k a s (n) 



K 2 7T 



(4.13) 



with the integration over the NP region of small k. It is simple to see that the regions 
non-collinear to the outgoing hard parton P 2 (away from (4.12)) lead to contributions 
which are even in b and then give subleading NP corrections. A correction linear in 
b is non-local in the conjugate variable Xi see [3]- 

To complete the calculation of the leading NP correction we use the standard 
procedure extensively discussed in [3, 16]. First we need to extend the meaning of 
the running coupling into the large distance region. This can be done by following 
the approach [5] in which one uses the dispersive representation 

s^ = r^»*sM (4.i4) 
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Moreover one needs to take into account the non fully inclusive character of the ob- 
servable H(x). This can be done in two steps, see [16]. One replaces the momentum 
k of the emitted soft gluon by k — > \/ n 2 + m 2 , with m the mass of the gluon decaying 
products. Then one multiplies the result by the Milan factor M.. Finally one obtains 

6B 12 (b) = -b • A NP , A NP = ~~~~ J dm5a cS (m 2 ) , (4.15) 

with 5a e s the large distance contribution of the effective coupling. See appendix D for 
the connection of A NP with the NP parameter already measured in 2-jet observables. 

Similar leading NP corrections are obtained for the B13 and B23 pieces. Since 
they come from the phase space regions collinear to the outgoing hard partons P2 
and P3, the complete NP correction from (4.8) is proportional to the total charge of 
the hard outgoing primary partons 

SB 5 (b) = - (C? } +Cf ) b ■ A NP . (4.16) 

In the present analysis we concentrate on the leading power correction discussed 
above. Large angle emissions and the ones collinear to the incoming parton generate 
next-to-leading NP corrections, proportional to b 2 ~ Q~ 2 . They appear not only 
in this B term, but also in the radiator 1Z (see (3.12)) and in the parton density 
functions V (see [5]). A complete study of higher orders NP corrections is still to be 
done. 

The final answer for J, including the leading NP correction in (4.16), thus is 
Mx; Q ab ) = 1°°^ ™s(b X P t ) F p {b) ■ {l - (Cf + Cf ) b ■ A NP | , (4.17) 

with J- p {b) the PT factor in (4.4). The final answer for the azimuthal correlation in 
(1.4) can be written in the form 

dx ax 

with b (and the outgoing hard parton charges) averaged over the PT distribution in 
(4.4). One can give an estimation of (b) by using the one- loop coupling approxima- 
tion. From appendix E we evaluate this mean at X — an d obtain 



1 MornV 4 2 + ^~C T 

(^=o~^r- Hr ' 7 = ^<7 r ln— 4— = 0.62 (n / = 3) . (4.19) 
Aqcd \ Q J Po 1 + 



Taking into account the two-loop coupling the exponent 7 is not a constant but 
smoothly dependent on Q, see [3]. The detailed numerical analysis will be done in 
the next section. 
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5. Numerical evaluation 

In the small x region, the PT distribution is obtained from (4.3), (4.4) with the PT 
part of the radiator in (4.1). At larger values of x one should use the exact PT result. 
In order to obtain a PT description valid for small and large x one needs to match 
the resummed and the fixed order distribution. This amounts in computing, to the 
due order, the coefficient function Cs{a s ) in (3.3). This is usually done by comparing 
the resummed expression with the exact result. 

We start from the integrated resummed distribution S PT (x) given by 




(bP tX )F P (b), (5.1) 



with JF p (6) the PT factor in (4.4) and a(y±) the Born cross section (see (3.17)). In the 
following we describe in detail the adopted procedure for the first order matching 
using a fixed order numerical program such as [17, 18]. Notice that, in order to 
control the a s scale, it is not sufficient to use first order matching but one needs 
to take into account also the two-loop analysis. To first order in a s the resummed 
distribution E PT (x) in (5.1) reads 

S PT (X) = \ (l + ^ (G 12 L 2 + G n (y±)L + cT(y ± )) + O (a s 2 )) , L = - ln X , 

(5.2) 

with G12 = —2Ct, Gn and c^ es given in (F.5) and a s = oc^(Q). Here the factor 
1/2 comes from the fact that the distribution (5.1) does not include the singular 
contribution proportional to S(x — 7r)- Only with this singular contribution the 
distribution is normalized to one. 

Following [3], we match the resummed expression in (5.1) with the exact first 
order result £i(x), which is here obtained using the numerical program DISENT of 
ref. [17]. We obtain 

^lt(x) = \ + (l + ^(ci - cD) (S pt (x) + SE( X )) , (5.3) 
with S PT (x) the resummed distribution in (5.1), the coefficient 

ci = lim (Ei(x) - G U L 2 - GuL) . (5.4) 
and the matching correction 

SZix) = I - G 12 L 2 - G n L - Cl )) . (5.5) 

Including the NP corrections we have: 

= \ + (l + ^(ci - cD) ((1 - A NP ((C 2 + C 3 ) • b)) S PT ( X ) + 6E(x)) , 

(5.6) 
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Figure 2: The azimuthal correlation distribution for s = 98400 GeV 2 , = 0.1, Q 2 = 
900 GeV 2 , y_ = 1.0 and y + = 2.5. The upper curve is the resummed PT distribution 
obtained from (5.1) (without matching), the lower curves are the PT and NP matched 
predictions obtained from (5.3) and (5.6) respectively. In the figure is also shown the first 
order prediction given by DISENT. 

where b is now averaged over the PT integrated distribution in (5.1). 

Here we report the result for the azimuthal correlation distribution with the 
following choices. We used the parton density function set MRST2001_1 [19] corre- 
sponding to a s (M z ) = 0.119. We also set a (2GeV) = 0.52 (see (D.23)), a value 
obtained from the analysis of two-jet event shapes [20]. 

In figure 2 we plot the full azimuthal distribution dH/dx obtained from (5.6) for 
xb = 0.1, Q 2 = 900 GeV 2 , y_ = 1 and y + = 2.5. The center of mass squared energy 
has been fixed at s — 98400 GeV 2 . In order to simplify the analysis (data are not 
yet available) we considered only the case of exchanged photon (see for instance the 
discussion in [7]). 

Together with this curve we plot the fixed order PT distribution obtained from 
DISENT and the PT resummed distributions obtained from (5.1) and (5.3). We 
notice that the PT resummed curve follows the fixed order one for not too small 
values of %■ Then, while the first order starts developing a singularity for \ — > 0, 
the resummed curve rises to a constant value. Actually the height of the 'plateau' is 
dramatically reduced after performing the first order matching. This is due to the 
large coefficient c x appearing in (5.3) (ci ~ —30 in the present case). 

The effect of the NP correction (5.6) is to further deplete the PT curve, in 
analogy with the EEC case [3]. 
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6. Discussion 

We presented the QCD analysis of E T EiC, the azimuthal correlation in DIS defined 
in (1.1)-(1.4). We considered small x, the back-to-back azimuthal angle, and finite 
t/2, the dijet resolution variable which selects dijet events and limits emitted jet 
transverse momenta not to exceed Q so that the exchanged photon is pointlike. This 
distribution is the generalization to DIS of EEC in e + e~. 

The QCD resummation is performed in the 1-dimensional impact parameter b 
conjugate to the soft recoil \P t = | ^ k X i\, with fc xi the out-of-plane momenta of the 
undetected secondary partons, see (2.7). We sum all < \n n+1 bQ and < \n n bQ terms 
in the exponent. The main complication with respect to the e + e~ case is that EtE-iC 
involves also initial state radiation which contributes both to the incoming parton 
evolution and to the observable H(x)- Due to QCD coherence the two contributions 
factorize, see (3.10). We have indeed that b^ 1 ~ x?t sets both the hard scale for 
the parton density functions V (instead of Q as in the cross section, see (3.17)) and 
the lower limit to frequencies in "radiation factor" A, a CIS quantity associated to 
the qqg hard vertex. See [7, 10, 12] for incoming radiation factorization in jet-shape 
observables in processes involving incoming hadrons. 

For the PT result at SL accuracy one needs to consider only the hard term H^ix) 
of the observable from the two outgoing partons in the qqg hard system, see (2.12). 
The general behaviour of dE PT (x)/dx (and its comparison with EEC) can be simply 
explained by observing that there are essentially two phase space configurations for 
the secondary partons: 

• all \k X i\ < xPt- This gives a contribution to gE pt \x) I dx which behaves as a 
Sudakov form factor (4.6) with a peak at small Xi see (1-6); 

• the other configurations, with some \k X i\ exceeding the soft recoil. This corre- 
sponds to impact parameter with all values. Their contributions to gE pt (x) / dx 
lead to a 1/ Y/a^-singular term (the integral of a Sudakov form factor) which is 
a decreasing function of x starting from a finite value at x — 0, see appendix C. 

The contributions from these two momentum configurations are comparable just at 
the Sudakov peak and for smaller x the 1/ A /a^-singular term is dominating. As a re- 
sult, dT l PT (x)/dx is continuously decreasing with x- Therefore the exponentiation of 
the one loop result in (1.5), with the characteristic Sudakov peak, is valid only above 
the Sudakov peak region. This is qualitatively different from the EEC behaviour in 
which the Sudakov peak is indeed manifest. 

We have assumed [5, 13] that the leading power corrections are obtained when the 
argument of the QCD coupling runs into large distance regions. Such a hypothesis has 
been successfully tested in 2-jet shape observables [13] and EEC [3]. The azimuthal 
correlation (as other 3-jet observables [6, 7, 10]) provides a further way to explore this 
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assumption. Indeed, the running coupling enters our result with the three different 
arguments K a b, the transverse momentum in the aft-dipole centre of mass frame. For 
each afr-dipole we evaluated the leading power corrections by taking n ab ~ A QCD , thus 
obtaining the result (4.16). The three soft regions K ab ~ A QCD are very different. 
Consider for instance the leading power correction from the 12-dipole which is coming 
from the momentum region (see (4.12)) 

k x = K x , \k y \ > \k v \ ~ \k x \ ~ Aqcd , (6.1) 

where k = k 12 while k t is the transverse momentum in the Breit frame. We have 
that NP region k ~ A QCD does not correspond to k t ~ Aq CD . One finds analogous 
differences comparing K i2 with K i3 or k 23 . This shows that analyzing EtEtC (as well 
as 3-jet observables) we can further check whether NP corrections are controlled by 
the soft arguments in running coupling. 

The strength of the leading power correction is controlled by the NP parameter 
A NP in (D.23) which is connected to the one in the thrust distribution, see (D.25). The 
leading power correction is similar to the ones in EEC: the expected 1/Q behaviour 
of the leading NP correction turns into 1/Q 1 with 7 a slowly varying function of Q. 
We have 7 ~ 0.62 for three quark flavours, see (4.19). 

In conclusion, this paper is a first step towards the full QCD analysis of azimuthal 
correlations in hard processes, which are among the most interesting and least studied 
observables in QCD. Resummation and NP corrections are both quite large. A 
detailed experimental analysis (yet missing) would therefore be important in order to 
test our knowledge of QCD dynamics and confinement effects in multi-jet ensembles. 
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A. Elementary partonic cross-sections 



For the elementary process (2.2) we may write the parton momenta Pi, P 2 as 
P 1 = ^(1,0, 0,-1), P 2 = Q(z ,0,T M ,z 3 ), 



(l-O(i-x) 



(A.l) 



a; a; 
in terms of the variables in (2.3). Distinguishing P 2 and P 3 according to (2.5), the 
variable £ is given in terms of Tm by 



1 



l - Wi - 



rT 2 

1 - X 



(A.2) 



and in terms of y 2 by 
1 — x 



2(2x- 1) 



i 4x(2a; - l)y 2 _ i 



0<x<^, 0<V2<^ 
3 1 n 2(1 -x) 2 

4 <JB<1 ' ° <m< xW^l) 



x 2 y 2 - (1 - a:) 2 
(2a;- l)(l-x + xy 2 ) ' 



3 1 2(1 -x) 
- < x < 1, 7- 

4 x(2x - 1) 



2 1-x 

<V2< 



X 

(A.3) 



The phase-space in terms of the variables (x, £) then becomes: 

1 



rr B < a; < , a;M 



i + y- 
l 



, 2/- < 1/3 



4y_ 



y- > 1/3 



(A.4) 



i 

2 ' 



X > 



4y+ 



We consider now the nature of the involved hard primary partons. We iden- 
tify the incoming parton of momentum Pi by the index r = q,q,g. Since (2.5) 
distinguishes P 2 from P 3 , in order to completely fix the configurations of the three 
primary partons, we need to give an additional index 8 — 1, 2, 3 identifying the gluon. 
Therefore the primary partons with momenta {Pi,P2,P3} are in the following five 
configurations 

{gqq} , or {gqq} 

{qgq} , 

{qgq} , (A.5) 
{qqg} , 
{qqg} ■ 



r = g, 5=1 
r = q,5 = 2 
T = q,5 
T = q,5 
T = q,5 
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Next we give the corresponding five elementary distributions da Tj sj, with / the 
fermion flavours. Since we consider only diagrams with the exchange of a single 
photon or Z , the individual partonic cross-sections may be decomposed according 
to: 



da 



T,8,f 



a 2 a s 



dxd^dQ 2 



[cf(Q 2 ) [(2-2y + y 2 )C r T 5 (x,0 + 2(l-y)C T L 5 (x,0 
+ Df(Q 2 )y(2-y)C T /(x^)}: 



(A.6) 



into transverse, longitudinal and parity-violating terms. Here y is defined by 

(Pq) Q 2 



y 



:pp p 



x B s 



(A.7) 



where P e is the momentum of the incident electron, and s is the centre-of-mass 
energy squared of the collision (we neglect the proton mass) . The flavour-dependent 
functions C^(Q 2 ) and D^(Q 2 ) show how the 7 and Z exchange diagrams combine: 



C f (Q 2 ) = ( 
Df(Q 2 ) 



Q 2 



2 2e f V f V e 

f sin 2 29 w \Q 2 + M 2 
2e f A q A 



+ 



(yf+A})(V?+Al 
sin 4 29 w 



Q 2 V 
Q 2 +M 2 J ' 



Q 2 \ W f A f V e A e f Q 2 V 
sin" 29 w \Q 2 + M 2 ) sin 4 29 w \Q 2 + M 2 ) 



(A.8) 



The functions C^ L / 3 (x,£) are the one loop QCD elementary square matrix elements 
for the hard vertex (there is no contribution of order O (a®) since we require trans- 
verse jets). We have [21] 



x 2 +e 



Cf = C F -4x£, 

x 2 +e 



r"i' 6 — c 



(l-x)(l-0 



+ 2(l + x0 



+ 2(x + 



(A.9) 

If we consider only photon exchange the index r is redundant, and the relevant 
elementary distributions dapj are 



, Cf = [x 2 + (1 - xf] 



Cf' 1 = 8x(l - x) , 



Of' 1 = , 



5,5 



s ^2 



dxdidQ 2 Q A f 

x 2 +e 



[{2-2y+y 2 ) C 5 T (x, + 2(1 -y) C 5 L (x, 0} , 

e + (l~0 2 (A.10) 



;i-^)(i-o 



+ 2(l+x£) 



, c^=[x 2 + (l-x) 2 } 



£(i-0 



C£ = CV • 4x£ , Ci = 8x(l-x) , C^ /L (x,0 = C^ /L (x, 1-0 . 
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B. Observable decomposition 

Here we decompose the observable H(x), defined in (1.1), in the case where we have 
two hard outgoing partons pi and p 3 accompanied by soft radiation fcj. Replacing the 
sum over hadrons in (1.1) by a sum over partons, and neglecting terms of O (k 2 /Q 2 ), 
we obtain 

H V0 = ~Q2 S (* - 7l )+ 2 -Q2- 6 (X ~ X23) 

fB 1) 

, n s^ Pt2k ti ; . / s Ptshi s( v 1 ' ; 

i i 

With Xa6 = 7T - |0 a6 |. 

In order to perform the calculation, we must introduce the definition of the event 
plane as the plane formed by the incoming proton direction np in the Breit frame 
and the unit vector n M which enters the definition of thrust major 

T M = max — y~]\ph • n M \ , n M ■ n P = . (B.2) 

h 

For dijet events with transverse momentum P t ~ Q we have T M = 0(1). Fixing 
the event plane as the yz plane, let the transverse momentum components of the 
outgoing partons be 

— * 

Pa = Ptfa, 1 - e 2 ) , Pa = -Ptfa, 1 - , hi = (k X i, k yi ) = fc ti (sin0 i , cos0;) , 

(B.3) 

with P t = \QTm- Conservation of momentum gives 

<p x = e 3 - ei = £ ^ , e 2 - e 4 = ^ ^ , (B.4) 

while the definition (B.2) of T M gives 

+ ^3 = - £ ^ - *(-**)] , e 2 + e4 = ^^. (B.5) 



Thus 



* = £^(-**) * = -£^(-**)- 

Expanding up to terms linear in the soft momenta, we get 

X23 = |e3 - ei| X2i = k + e i - 0i| X3i=|0i-e3| 
P<2 = - e 2 ) Pt3 = PtO- - e 4 ) 



(B.6) 



(B.7) 
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Hence the observable H(x) is given by 
IP 2 

H{x) ~(l - £2 - Cd) (X - |es " cil) + % - *)] 

+ 5^ fc « - k + e i - &|) + <Kx - \<t>i ~ e 3 |)] • 
^ i 

The first term is generated by the two hard partons, but has contributions both from 
the underlying hard momenta and from the soft parton recoil. The second term is 
the direct contribution from soft radiation. 

Let us therefore take the recoil contribution from the first term (that proportional 
to e 2 + e 4 ) and add it to the soft term: 

op2 9P r 

H (x) = -7W " I 63 - ei ') + ~ *)! + ni E k A 5 ^ ~ \* + 61 " 

+ <Kx - l& - e 3|) - I cosfa\5(x - |e 3 - ei|) - | cosfa\5(x ~ n) 

Note the cancellation of collinear singularities: when ki is collinear to p\ we have 
fa = ei, which generates cancellation between the second and third delta- functions 
in the soft term above, and between the first and fourth. (We have cos fa = 1 to this 
accuracy.) Similarly, when ki is collinear to p 2 we have 0« = n + which generates 
cancellation between the first and third delta-functions, and between the second and 
fourth. (Now we have cos0j = — 1.) 

The above is valid for all x U P to x — tt, but we are only interested in small x 
for the resummation and leading power correction, with larger values of x given by 
matching with the fixed order calculation. For small x, we may write 

IP 2 2P r - 

H(x) * -^-6( X -IM+qiY, hi [6(x ~ \4>i ~ M - I cosfa\5( X - \fa 

(B.10) 

where instead of fa we now use fa, the azimuthal angle between ki and either the 
hard parton p 2 or p 3 : 

<Pi2 = fa - ei for h near p 2 , ^ ^ ^ 

fai = it - fa + e 3 for ki near p 3 . 

The cancellation of the collinear singularity is still manifest as we take fa — > 0. 
C. DL approximation 

Here we give an instructive simplified study of the behaviour of the small and large-6 
contributions to the distribution We use the DL radiation factor evaluated at 
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Pt — \Qi with a fixed value of the running coupling (a s = 0.119) and without parton 



2 

distribution functions: 



^dl(6) = e- R ^ , J R DL (6) = ^Iln 2 (6P i ), b = be™. (C.l) 



7T 



We separate the 6-integral into a soft and hard contribution as indicated in equation 
(4.5), and discuss each contribution in turn. 

C.l Large b (soft) contribution 

The large b or soft contribution is that arising from multiple gluon emissions with all 
\k x \ < xPt, whose cumulative effect is a hard parton recoil of xPt- Thus we expect 
to see in this contribution the usual Sudakov behaviour. We obtain 



J(+) { X ) = 4- ^7 MPtbx) e- RDL(b) 
dx Ji Kb 



m 
dx 



where the function R' is 



e-v** sec f R' lg (-1)" e -(^+ih E 



2T(1 + R>) ir ^ (2n+l)\2n + l - R' 

\ J n _Q \ / 



(C.2) 



J?=^hi, (C.3) 

X 

and the quantity in brackets is a SL function which equals \ at DL level. The result 
is well-defined and convergent for all R ', since the poles from the sec function are 
cancelled explicitly by the second term in the brackets. (If however we attempt the 
same procedure on the integral over the whole b range, we obtain only the first term, 
which diverges at R' — 1.) 

The difference between the approximation and the exact value of the integral is 
formally beyond SL and is compensated for by the matching with fixed order. The 
Sudakov behaviour is clearly seen, along with the peak at R' ~ 1 (see figure 3). 

C.2 Small b (hard) contribution 

The small b or hard piece is that arising from two or more hard gluon emissions, with 
\k x \ > xPt, but which combine with opposite signs to give the hard parton recoil of 
only xPt- We can expand in powers of X'- 

l ( -\x) = j- [^-sm(P t bx)e- R ^ 

dx Jo ^ , . 

nd X ^ (2n + l)! J 
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Figure 3: Hard, soft and total contributions to T{x) in the DL approximation for the 
azimuthal correlation. 

and then integrate to give 

1 d j2 (-l) n (xe" 7B ) 2n+1 e^ 1 ^f^^l 2n + l 

where the quantity a and the function are 

2ouG 



z<-> M =~y, ( - i) ; 2 (xe :;T" e ^h* (c ,> 

(2n+l)! Va V. X V® J 



7Y 



= / dte^* 2 . (C.6) 



The series is rapidly convergent since as n increases the argument of the $ function 
becomes more negative. For this contribution there is no Sudakov behaviour, but 
rather at x = the distribution increases to a constant that behaves as 1/ y/a^. 

In figure 3 are plotted the two contributions to I as well as the total. Notice 
that at very small x (i- e - with R' > 1) the soft radiation is suppressed and the 
distribution is due to hard emission. Soft radiation contributes around and above 
the Sudakov peak at R' — 1. 

C.3 Comparison with EEC 

The case of energy-energy correlation in e + e~ is similar to the above, except that 
instead of Ct the total colour charge of the hard partons is 2Cp, and the integration 
measure is substituted according to 

* sin(P t 6 X ) — > Q X db J, (Qbxj . (C.7) 
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Figure 4: Hard, soft and total contributions to I(x) in the DL approximation for the 
EEC distribution. 



In particular the 2- dimensional phase-space ensures that the distribution peaks and 
falls linearly to zero, while the fact that the SL approximation is valid up to R' = 2 
means that it is a good approximation for in the region of the peak, up to matching 
with fixed order. 

In figure 4 are plotted large and small-b contributions and the total for this 
distribution. The large-b contribution dominates until well below the peak. This 
very much simplified study shows the qualitative behaviour of the distribution in 
6-space, and gives a helpful insight into the physics of the azimuthal correlation. 



D. Soft contribution 

In this appendix we study the perturbative and non-perturbative contributions to 
the quantity B$(b) given in (3.13). We consider each dipole separately. 

D.l Dipoles 12 and 13 

Consider 

/H 3 k hi 
— Wl2 (k) -| [cos(W) - I cos0|] . (D.l) 
mo r t 

Denoting by pl,p 2 and k* the momenta in this system, we introduce the Sudakov 
decomposition 

P5 = %1, 0,0,1), pj = %(l, 0,0,-1), k* = ap* 2 + Ppl + K, (D.2) 
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where the scale Qf 2 = 2{pip2). Here the two-dimensional vector k = (k x , k v ) is the 
transverse momentum orthogonal to the 12-dipole momenta (k 2 = k 2 2t ). We have 
then 

,,. aJn 2 ) d 3 k d 2 nda k 2 n . 

W12 (fc) = 5~ , = , a > —J- . D. 3 

We must also express the Breit-frame components of in terms of these variables. 
The Sudakov variables a and (5 are Lorentz invariant, while the two-component vector 
k transforms into a 4- vector with components 

Breit _ P t Breit _ Breit _ Breit _ _EL K (T) A\ 

Thus the transverse momentum components of in the Breit frame are 

k x = K x > = aP t + K y • (-D-5) 

We take as the azimuthal angle between p 2 and k — this neglects the contribution 
from k near to y>3, which contributes beyond SL and so is included by matching with 
the fixed order result. Therefore 

- k v k v aP t , aP t , 
cot = y- = — H = cot H CSC , (D.6) 

•^3? ^ 

and the 12-dipole contribution is 

[ dn 2 d<p da a s (n) /tsin0 f - -, 

Bi 2 (6) = / — - . T cos(W) - | cos 0| . D.7 

J K 27T Ct 7T it sin 

Changing integration variable from a to gives 



n , f*d<f> . 2A [+" rf0 cos(P t 60)-|cos0| 

B\2\o) — —fr / dK,a s {K) / — sin -* J 



nP t J ' 7 t J4> m sin(0-0) sin 

COt m = COt H CSC , COt (j)M = COt + — 2~ CSC , 



(D.8) 



which is clearly beyond SL, since both the soft logarithm (k — > 0) and the collinear 
logarithm (0 — > 0) are cancelled. The 13-dipole is treated similarly. 

D.2 Dipole 23 

The 23-dipole is analysed along the same lines as the 12, except that there are some 
minor differences. Again we introduce the Sudakov decomposition 

^ = %(1, 0,0,1), p* = %(l, 0,0,-1), k*=ap* 2 + (3 P ; + K. (D.9) 
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We must again express the Breit-frame components of k in terms of these variables, 
but now the two-component vector k transforms into a 4- vector with components 



Breit _ Breit _ Breit _ /-. _nt\ K ..Breit _ ' 1 . 

~ Q(l—x) ~ y ~ S'^f > K z ~ Q(\—x) y ' 

(D.10) 

The transverse momentum components of k in the Breit frame are then 

k x = K x , k y = (a-/3)P t + (l-2£)K y . (D.ll) 

We take as the azimuthal angle between p 2 and k in the region where a > j3, 
and as the angle between p 3 and k when a < (3. The neglected contributions are 
beyond SL and so are included by matching with the fixed order result. We take 
each hemisphere separately: in the right hemisphere (a > (3) 

cot0 = (1 - 2g) cot0 + ^ ~ ® P * esc , (D.12) 

K 

and the 23-dipole right hemisphere contribution is 

„2 

l a>j3 k 2 2-n a n P t sin i 
Changing integration variable from a to gives 

d(j) cos(P t 60) — | cos i 



d / 7 x f dn 2 d6 da oiJk) Ksin0 r . „ , -. 

- I -^i-^r • ^ t cos(w) " 1 cos ^ • (D - 13) 



V^0) sin 2 



- 2 P 2 

f((f>,<f>) — [sin cos — (1 — 2£) cos sin 0] + 4— |- sin 2 

^23 

cot m = (1 — 20 cot H CSC , cot 0M = (1 — 20 cot , 

(D.14) 

which is clearly beyond SL, since both the soft logarithm (k — > 0) and the collinear 
logarithm (0 — > 0) are cancelled. The left hemisphere is treated similarly. 

D.3 NP contribution 

We calculate the NP contribution using the standard procedure, see [5]. 
Applying this to equation (D.8) gives the NP contribution to B\ 2 as 

,u tu\ M [a 2; / 2^ ~d f°° d ^ [* d< t> ■ 2, 

dB 12 (b) = — / dm ba efi {m )-— / = / — sin 

7T 7 dm 2 Jo PV/t 2 + m 2 Jo ^ 

d(f) cos(P60) - | cos 0| 



X 



sin (0 — 0) sin 2 



t , Pt 7 Pv^ 2 + m 2 

COt m = COt H ^=:^= CSC , COt 0M = COt H -75 CSC 

AHm 2 Q12 



(D.15) 
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Here Ai is the Milan factor [16] which accounts for the non fully inclusive nature of 
the observable. 

Doing the m-derivative and ^-integral gives 

2M [a 2 , fd^ . 2 . [* M d<p cos(W) - | cos 0| 



s sin(0 — 0) sin 2 



dB 12 {b) =— — / dmda cS (m ) — sin / 
TT-n J Jo 11 J4, 

7 i i 7 PfTn 

cot m = cot 0-1 csc , cot <pM = cot + » CSC . 

m gf 2 



(D.16) 



We wish to find the leading term in the expansion of this integral in the range 
m <C b^ 1 <C P t . Writing = u/(bP t ) and expanding the integrand yields 

.„ /L , 2M f, . , 2^ ■ i frnj COSM-1 

o B 12(b) = / amoa c g{m ) / — sm0 / bP t du - h • • • 

kPi J Jo ^ Jo u 2 , s 

= —b ■ A NP + • • • A NP = ~~~ I dm5a cS (m 2 ) , 

where the dots involves terms which are of higher order in b. The 13-dipole gives an 
identical result. 

The same procedure applied to equation (D.14) gives the NP contribution to B^ 

as 



6B*{b) J± J dm 2 5a eS (m 2 )^- 2 J 



dn 2 



+ m 2 Jo vr 



' • 2 
■ sin . 



x rf0 ^ cos(P f &0)-|cos0| 

p _ 
COt m =(1 — 2£) COt H -^=^= CSC , COt (f>M = (1 — 2£) cot . 

V k 2 + m 2 

Doing the m-derivative and ^-integral gives 

2M f a , , 2 , . 2i f* M d(j) cos(P t 60) - I cos 0| 



SB^(b) = 2 ^L f dm5a cfi (m 2 ) [ — sin 2 / 

TT-Pt 7 Jo 11 Jd, 



4> m v7(070) 

p 

COt m =(1 — 2£) COt H CSC , COt 0M = (1 — 2£) COt 

m 



sin 2 



(D.19) 



We wish to find the leading term in the expansion of this integral in the range 
m <C b^ 1 <C P t . Writing = u/{bP t ) and expanding the integrand yields 

COSM — 1 

ir (D.20) 



SB^b) — — — / dm 8a e s(m 2 ) / — sin0 / bP t du 
•• ' / ./ .Ai Jo 



-b • A NP + 
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where, as before, the terms in the dots are of order b 2 . The left hemisphere gives an 
identical result, and so the total is 



SB 23 (b) = -2b- A NP + --- (D.21) 

Combining the three dipoles then gives the leading non-perturbative contribution 

5B s (b) = - (cf + Cf ] ) b ■ A NP . (D.22) 

The complete treatment of this correction requires the definition of the coupling 
in the infra-red region (this will be done by the dispersive method [5]), the analysis 
of contributions which are non fully inclusive, leading to the Milan factor [16], and of 
the merging PT and NP contributions to the observable in a renormalon free manner. 
This can be done in the usual way and the NP parameter is given by 



A NP = M^m {aoM ~ »s - ( ^ — + ^ 4 1 



7T Z 



27T \ Pi A) 



(D.23) 



where 



3 - UNc 

PO — n 



2n f 



(D.24) 



The K factor accounts for the mismatch between the ms and the physical scheme 
[15] and ao(/ i /) is the integral of the running coupling over the infra-red region, see 
[5]. We have that, at two loops, A NP is independent of ///. Assuming universality 
of the coupling [5, 20] this NP parameter is the same as already measured in 2-jet 
observables. For instance, the shift for the r = 1— T distribution is 



%{t) = ^(r-A T ) , A r = Oc r A NP , c T 



(D.25) 



where Cf enters due to the fact the 2-jet system is made of a quark-antiquark pair. 



E. Estimates of (b) using 1-loop coupling 

We consider using the 1-loop running coupling 

a s (k) 



A) ln(fc/A) 

The DL function r(b,Q) defined in (4.1) then becomes 



r(b,Q) 



A) 



(E.1) 



b< 



A ' 



(E.2) 
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and we impose r(b,Q) — oo for b > |-. We also use parton density functions deter- 
mined by the 1-loop evolution equation 

dV N (k) a s (k) 



dink 7r 



-lN-V N {k) (E.3) 



with 7n the leading-order anomalous dimension matrix, and VN(k) the moments of 
the parton density functions at scale k. This equation has solution 

We define the integrals 

roc 

ljfl(N)= dbb n V N (2b- 1 )e- Rs{l) , b = be~< E , (E.5) 
Jo 

using the radiator given in (4.1). The choice of 2b~ l rather than b^ 1 as the fac- 
torisation scale is merely for convenience — the difference is beyond SL. Evaluation 
yields 

li n ) ( N )=(-^) I " I i f 4^ 1 (E-6) 

with 



—) Sr— 

The average value of 6 is therefore, at fixed N, 



x 



(, 4 ^(«) , x 4 r (S) . d S) Q ( a S) 



(E.7) 



(i> " 4^ ~ — UJ (E ' 8) 

with a coefficient a«s that depends on N. The exponent is ^-independent and then 
we obtain the DL estimate in (4.19). 

F. Formulae for numerical analysis 

In this appendix we report some analytical formulae needed to perform the numerical 
analysis described in section 5. 
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Within SL accuracy the PT radiator (4.1) can be written as 

n (S)AS) 

Rs(b) = C T ri (5) - C T In 2 r'{b) + £ Cf In r 2 (b) , 5 = be" , (F.l) 

a 

where , ( { a 6) and are introduced in (4.2), while ri,r 2 and r' are given by: 

n(6) 



f Q dk 2a s (k) Q 8tt ln , ^ . . 4K / . „. £ \ 



4/3, A 2 ln(l-<)+^ «... 

r2 ( 6)= = * N1 _ 0i r , (5) = ML!) ln5Q - 4 



(F.2) 

In the above expressions /3o can be found in (D.24) and (5\ is given by: 

A = - 5CAHf - 3CFHf . (F.3) 

3 

The coupling a s is in the ms scheme and the constant K relating the physical 
scheme [15] to the ms is defined in (D.24). 

Actually (F.2) makes sense only for I < 1, so that, as usual (see [3]), we impose 

e -m = o ; for Iq > exp ( ) . (F.4) 

We report also the expression for the coefficients Gn and c^ cs introduced in (5.2): 
, , /*+ / da p \ 4P2 , , 



^^°» 2P, + Q, dlnQ j ' 

(y±) = a (y±)J2 J ~l_ «(^) q^(t'«) 



In^Vc^ln^ I 27 ^ ln M 



-C T ^-2C T ln 2 § 
6 Q 



(F.5) 
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